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Abstract 

An experimental scheme is suggested that permits a direct measure of en- 
tanglement of two-qubit cavity system. It is articulated on the cavity-QED 
technology utilizing atoms as flying qubits. With this scheme we generate 
two different measures of entanglement namely logarithmic negativity and 
concurrence. The phenomenon of sudden death entanglement (ESD) in a 
bipartite system subjected to dissipative environment will be examined. 
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1. Introduction 

Entanglement is an essential feature of quantum mechanics that permits 
basic peculiarity between classical and quantum physics. It is at the heart of 
many applications of quantum information science, including quantum tele- 
portation [Jfl] , quantum dense coding [5j], quantum cryptography H, and 
quantum computing Entanglement can reveal the nature of a nonlocal 
correlation between quantum systems that have no classical interpretation. 
Recently, the use of the cavity QED in quantum information processing be- 
comes more interesting [8-12], and the entanglement generation and nonlo- 



cality test of two cavity fields have been explored [l3Hl7||. Real quantum 
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system are unavoidably subjected to their environments, and these recipro- 
cal interactions often result in the dissipative evolution of quantum coherence 
and loss of useful entanglement. Decoherence may guide to both local and 
global dynamics, which may incite the eventual deterioration of entangle- 
151 ] . Particularly, Yu and Eberly have shown that the entanglement of 



ment 



bipartite systems can decay to zero abruptly, in a finite time which depends 
upon the initial preparation of the atoms, a phenomenon termed entangle- 
ment sudden death (ESD)]18j and was recently observed in two sophistically 
designed experiments with photonic qubits 3]and atomic band [io|. Fur- 
thermore, it has also been observed in cavity QED and trapped ion systems 



2l| . On the other hand, the phenomenon ESD has provoked many the- 
oretical investigations in other bipartite systems involving pairs of atomic, 

multipartite systems [26] and spin chains 



ohotonic, and spin qubits |22H25 
" 2 

la 



27 
16 




In addition, ESD has also been explored for different environments 
However, numerous investigations on ESD in a variety of 
systems have been done so far, the question of ESD in interacting qubits 
remains yet open 32 . 

On the other hand, several methods have been proposed for detect- 

ing and measuring entanglement without a full reconstruction of the state. 
These methods, although much simpler than the full state reconstruction, are 
not completely free of experimental difficulties, as they require either con- 
trolled unitary operations or some prior knowledge about the quantum state 
in question, or they can detect entanglement but not measure its amount. 
Furthermore, for two qubits the defining measure of entanglement is concur- 
It is a good measure of entanglement in every sense and direct 



rence 
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measure of the concurrence of two-photon pur e entangled state was confirmed 
experimentally using linear optical means [39[ . Nevertheless, it exists another 
computable measure of the entanglement called negativity 4(3| , and thereby 
fill an important gap in the study of entanglement. It can be regarded as 



a quantitative version of Peres' criterion for separability [41]. For a mixed 



quantum states, the two measures are different. We will show some analytic 
relations between the two previous measures of entanglement for the pro- 
posed system. 

The aim of this paper is two-fold. First, we propose an efficient scheme for 
quantum teleportation to generate entangled number states of two-bipartite 
system under the influence of dissipative environments. Second, we present 
two different computable measures of entanglement namely, the logarithmic 
negativity and the concurrence and we compare their amounts and agree- 
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ments for different reservoirs. Thus we investigate the problem of ESD for a 
this proposed scheme. 

The format of this paper is as follows. In section [2] we investigate the model 
for two-bipartite system with a simple dissipative reservoir and formulate 
their dynamical evolution by solving the master equation of motion. In sec- 
tion [3] we present the theory of two different measure of entanglement of 
two-bipartite system subjected to dissipative environments (the logarithmic 
negativity and the concurrence). In section H] we study the entanglement 
dynamics of two-bipartite system in vacuum and thermal reservoirs. Then 
we compare the amounts of the two computable measures of entanglement 
for different initial entangled states. We find that for thermal reservoirs the 
ESD always exists but for the vacuum reservoirs the ESD can be exhibit 
with some entangled states. Finally in section [5] we conclude with a general 
remarks and future outlook. 



2. Model 



Recently, Zubairy et all [43( have proposed a new scheme in their inves- 
tigation of the quantum disentanglement eraser. In this simple scheme , the 
concurrence can be directly measured from the visibility for an explicit class 
of entangle states. We propose here the same scheme but with some modi- 
fications. A two-level atom with the upper level |e) and the lower level \g) 
passes consecutively through cavity A, a field damping region and a cavity 
B as shown in figure IA.1I The atom is initially prepared in the excited state 
|e) and the decay of the radiation field inside a cavity may be described by 
a model in which the mode of the field of interest is coupled to a whole set 
of reservoir modes. We assume that initially the two cavities are in vacuum 
state |0) and the atom always leaves the setup in the ground state \g). 
In the interaction picture and the rotating-wave approximation, the Hamil- 
tonian is simply 

j 

j 

where cla(b) and a A , B \ are the destruction and creation operators of the mode 
of the electromagnetic field of frequency v. b^ B ^ and b^ A ^ are the modes 
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of cavity A(B) of frequency Uj which damp the field and is the coupling 
constant of the interaction between the electromagnetic field and the cavity. 
From the general analysis of system-reservoir interactions, with the Hamilto- 
nian (jTJ), we can obtain directly the master equation for the reduced density 
matrix for the filed in the cavities as 42 



P(t) 



„ 7 « _ , 



i=A,B 



a A 



i=A,B 



7 



(0 



-m 



aia\p{t) - 2a\p(t)a,i + p(t)aja| 



(2) 



where 7« is the decay rate in the cavity, and rii(i = A,B) are the average 
number of quanta at frequency v in the thermal reservoir which surrounds 
the cavities A and B. If the reservoirs are at zero temperature, n, = 0, and 
the remaining terms are due to vacuum fluctuations. 

In the general case, we consider the field states in Fock basis in two identical 
high-Q cavities A and B that represent a bipartite system containing the 
entangle field as 

|*)ab(0) = ai|0) A |0) fl + a 2 |0) A |l)s + a 3 \l) A \0) B + a 4 |l)A|l)fl (3) 



where aj(i = 1, 2, 3, 4) are the probability amplitudes with J2t=i l a «| 2 = 1- We 
use the basis ( |1> = |0) A |0) B , |2) = \0} A \1) B , |3) = \1) A \0) B , |4) = |1)a|1)b) 
to define the density matrix of the two qubit system. The equations of motion 
in terms density matrix elements can be obtained using the master equation 
[2] and with their solutions in the general case are given in the Appendix A. 



3. Degree of entanglement 

To study the effect of interaction among the two-bipartite on decoherence 
we have to investigate the dynamics of two-bipartite entanglement. In order 
to compare the degree of the entanglement restrained to quantum state, we 
will use two entanglement measures , i.e., logarithmic negativity and concur- 
rence, to describe the degree of entanglement for any bipartite system. Both 
measures satisfy necessary conditions for being good measures of entangle- 



ment. The logarithmic negativity [40l. |44j for two-bipartite system is defined 
by 

A/" = log 2 ||p Ts ||i (4) 
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where p Ts is the partial transpose of a state p in d <8> d! (d < d') quantum 
system and ||.||i is the trace norm that can be read as 

Hp^lK = 1 + 21^1 (5) 

i 

where /ij are the negative eigenvalues of p Tj3 . For pure states, Af = for 

unentangled states and M = for the maximally entangled state. 

We will also consider another important measure of entanglement that is the 



concurrence 



38, 45 



C(t) = max{0, V% - V A 2 - V A 3 - V^) (6) 

where A's are the eigenvalues of the non-hermitian matrix p(t)p(t) arranged 
in decreasing order of the magnitude. The matrix p(t) is the density matrix 
for the two-bipartite and the matrix p(t) is given by 

pW = «®<)p*W«®<) (7) 

where p(t)* is the complex conjugation of p(t) and a y is the Pauli matrix 
given in quantum mechanics. The concurrence fluctuated between C = 
for a separable state and C — 1 for a maximally entangled state. The two 
measures of entanglement are different for mixed quantum states. 
Here we will consider some interesting initial entangled states for the two- 
bipartite which can be prepared and have potential applications in the quan- 
tum information processing tasks [3-0]. 

We will start by the investigation of the EPR-states which are concepts in 
quantum information science, a crucial part of quantum teleportation and 
represent the simplest possible examples of entanglement. 

1. Assume that the initially entangled state of the field in two cavities to 
be in a NOON state given by 

|*)ab(0) = a a |0) A |l) B + a 3 |l)A|0) B (8) 

This kind of state can be generated by passing a two-level atom initially 
in the excited state through the two empty high-Q cavities. The inter- 
action times of an atom with two cavities are chosen to be such that we 
have a 7r/2 pulse in the first cavity and a ir pulse in the second cavity 



461 ] . The initial logarithmic negativity and concurrence are given by 
JV(0) = log 2 (l + 2|p 23 (0)|) = log 2 (l + 2|a 2 a 3 |) (9) 
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C(0) = 2|p 23 (0)| = 2|a 2 a 3 | (10) 

where the two quantities are related for this case by Af(0) = log 2 (l + 
C(0)). It is evident, the solution of this equation gives A/"(0) = C(0) = 
1 which corresponds to the case |a 2 | = | «s [ = In (Figure IA.2j) 
the variation of the initial values of the logarithmic negativity and 
concurrence in terms of |<Z2 1 is presented where 1 0,2 1 2 + 1 0,3 1 2 = 1. 
Using the solutions of Appendix A, it can be shown that the density 
matrix p(t) can be read as 



P(t) 



{ Pu(t) \ 

p 22 (t) f>23(t) 

P* 23 (t) Pss(t) 

\ Pu (t) j 



Then, we can calculate the negativity defined by ffT4|) for the two- 
bipartite. We find 

M(t) = max[0, 



log 2 



1 - p n (t) - p u (t) + v/[p n (t) -p 4 4(t)] 2 + 4|p 23 (t)| 2 



while the concurrence defined by ([6]) is given by 

C(t) = max (0,2 \\p2s{t)\ - y/pn(t)p u (t) 



(12) 



(13) 



As they decay, they get entangled with the environment, slowly losing 
their coherence and purity over time. 

For the case of vacuum reservoir, we can see from the solutions of the 
Appendix A that pu{t) = 0. Then the relation between concurrence 
and logarithmic negativity can be given by 



Af(t) = max (0,log 2 [l - pn(t) + vViW 2 + C 2 ] ) (14) 
C = 2|p 23 (t)| (15) 

where we can not observe in this case the ESD for any initial states. 
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2. Consider now the initially entangled two-bipartite to be in a another 
EPR-state given by 



|^>AB(0) = a 1 |0) A |0) B + a 4 |l>A|l>B 



(16) 



This kind of states can be prepared by swapping the vacuum and one- 
photon state in the cavity A of the state (jSJ) discussed above [2lj, |47 . 
States like this have been realized in experiments with trapped ions 



48l |49| . The initial logarithmic negativity and concurrence read as 

JV(0) = log 2 (l + 2|p 14 (0)|) = log 2 (l + 2\ ai a 4 \) (17) 

C(0) = 2|p 14 (0)| = 2|o 1 a 4 | (18) 
where the maximum value of A/"(0) and C(0) are equal 1 for |ai| 



where \a\\ A + |a 4 | 



1. The initial values of logarithmic 
negativity and concurrence have the same behavior (Figure IA.2[) in 
terms of |ai|. 

Exploiting the solutions of the Appendix A, the density matrix p(t) 
can have the form 



P(t) 



( Pn{t) p 14 (t) \ 

p 22 (t) 

o p 33 (t) 

V pIM o o PM (t) j 



(19) 



We can show that the logarithmic negativity defined by (JSJand the 
concurrence defined by dSJ) for the two-bipartite in this case are given 
by 



M (t) 



maxl 



log 2 



1 " P32(t) ~ P33(t) + V / [P2 2 (t)-p33(t)] 2 +4|pi 4 (t)| i 



and 



C(t) = max 0,2 |p 14 (t)| - y'ftHj 



(20) 
(21) 
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As is evident form the solutions in Appendix A, in this case p 2 2(t) — 
p 33 (t). Then the relation between concurrence and logarithmic nega- 
tivity can be given by 

Af(t) = max (o,\og 2 [l + C(t)]J (22) 

C{t) = max(0,C(t)) (23) 

where 

C(t) = 2[\p u (t)\-p 22 (t)} (24) 

which makes clear that the logarithmic negativity is greater then the 
concurrence, when C decreases in terms of t, except for the initial value 
which is one for both of them for |ai| = | a-4 1 = and the value for 
t — > 00 which is also zero. 



4. Entanglement dynamics 

We will explore the time evolution of entanglement of the previous en- 
tangled states of two-bipartite system exposed to either vacuum or thermal 
reservoirs. In figure IA.3I we plot the logarithmic negativity and the concur- 
rence for the first EPR state discussed previously in (i), for different initial 
states for vacuum reservoirs. Note that for vacuum reservoirs n = 0, i.e. 
at zero temperature , no ESD is observed. A/"(t) and C(t) monotonically 
decrease to zero as t — > 00. On the other hand, the two measures give differ- 
ent values for entanglement in vacuum reservoirs. Generally, the logarithmic 
negativity takes smaller values than the concurrence 150], except for initial 



value for some initial states figure (IA.3^ b) . this can be seen from the equa- 
tions dm ITS]) where the concurrence is proportional to the population p 2 3(t) 
which monotonically decreases, while the logarithmic negativity depends on 
the population pn(t) which growths rapidly to rich a maximum value but 
the population p44(t) remains zero. In figure IAT41 we present the two mea- 
sures of entanglement for the thermal reservoirs in the case \a 2 \ = 1 0,3 1 = 
and for different values of n. Here the ESD is observed and as the tempera- 
ture increases, the sudden death time decreases. On the other hand, we can 
observe the perfect correspondence between the logarithmic negativity and 
the concurrence for the thermal reservoirs while in the vacuum reservoirs the 
concurrence exceed the logarithmic negativity which decays faster than the 
concurrence where The population pa{i) starts to manifest which restrains 
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the increasing population p n . The sudden death time is identical with the 
two measures of entanglement. Indeed, for the 2 two-dimensional systems, 
logarithmic negativity and concurrence are real entanglement measures: they 
do not increase local operations and classical communication, and vanish if 
only if the state is separable. 

In figure IA.5I we plot the dynamical evolution of the entanglement when the 
second EPR state (ii) is considered, for | and for different val- 

ues of n. For this case no ESD in the vacuum reservoirs, but when ai < 04 
we can see a finite-time disentanglement, which means that the major con- 
tribution of this state is responsible for the ESD. Thus, we can conclude 
that locally equivalent pure states with the same entanglement perform very 
differently during the time evolution and simple local unitary operation acts 
on the initial state can give naissance to ESD. As the mean photon numbers 
in cavities increase, the ESD is always observed and the finite-time disen- 
tanglement persists as shown in figure IA.6I Contrary to the first EPR state, 
the logarithmic negativity exceeds the concurrence and they are in good co- 
incidence. Furthermore, the sudden death time is the same with the two 
measures of entanglement. 

On the other hand, it is obvious that the two measures of entanglement are 
very closed. The sudden death time is the same in the two measures of en- 
tanglement, and the logarithmic negativity predicts the same behavior of the 
entanglement as the concurrence. When the average thermal photon number 
is different from zero, we can see that the ESD always occurs whatever the 
initial entangled states are and no matter how the nonzero average thermal 
photon number is. Furthermore, the two proposed measures of entanglement 
coincide and give a good prediction of the sudden death entanglement of the 



two-bipartite. This is consistent with the findings in [15|, [16 



5. Conclusion 

The delicate aspect of entanglement has been subjected to many quantita- 
tive studies and has guided to interesting results. Understanding the physical 
meaning of entanglement measures continues, however, a major defy. In this 
work we have shown that the scheme we propose here allows direct measures 
of entanglement: logarithmic negativity and concurrence of a two-bipartite 
cavity system. The sudden death time depends on the number of photons 
in the cavities and the temperature of the reservoirs. It increases with in- 
creasing the number of photons in the cavities and decreases with increasing 
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temperature of the reservoirs. The proposed scheme only involves system- 
reservoir interactions corresponding to the decay of the radiation field inside 
a cavity (field damping). This operation is based on the elementary exchange 
of energy between system and reservoir that is thus assumed to consist of 
the simultaneous creation of a quantum excitation of the system with an- 
nihilation of a quantum in one mode of the reservoir, or reverse process. 



These operations have been demonstrated experimentally |51[ and therefore 
our proposed scheme can be realized within the present cavity-QED tech- 
nologies. As a future perspective il would be interesting to study the effect 
of other type of environment like a squeezed reservoirs. Further, It would be 
motivating to enlarge our point of view to the two-bipartite entanglement of 
high dimension for the study disentanglement and loss of decoherence, where 
the Wootter's criterion is not applicable and we can use logarithmic negativ- 
ity which is necessary and sufficient condition for 2X3 and 3X2 systems. 

Appendix A. Equations of motion of the density matrix elements 
and their solutions for vacuum reservoir 

The equation of motion of density matrix elements for the general state 
(Eq. [3]) are given by 

pn(t) = i-n A l {A) ~ n Bl ^) Pll {t) + (n A +) 7 (A) P33 (f) + (n B + 1) T (B W*) 
P12W = -i\(2n B + 1) 7 (B) + n Al ^] Pl2 (t) + (n A + 1) 7 (A W(*) 

Mt) = -[^(2714 + 1)7^+ W B Vl3W + (^B + 1)7 (B W*) 

Mt) = -[(nA + \h {A) + l {B) (n B + 1 -)]p u (t) 

P2i(0 = -[(^(2nB + l))7 (B) +nA7 (A) ]p2iW + (^ + l)7 {A) P43W 

p 22 (t) = -[n A l {A) + {n B + l)7 (B) ]P2 2 (t) + in A + l) 7 (A) p4 4 (t) + W 

P23W = -[(^ + ^)7 (A) + (^ + ^)7 (B) ]p23(t) 

hS) = ~[(n B + 1))7 (B) + (n A + ^) 7 {A) ]p24(t) + n B j^p 13 (t) 
P3i(t) = -^(2rr A + l))7 (A) +W B) ]P3iW + (^ + l)7 (B) P4 2 W 
P32W = '[(n A + \h {A) + (n B + ^]p 32 (t) 
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P2a(t) = -[(n A + 1)7 (A) + n Bl {B) ]p 33 {t) +n A 7 {A ViiW + (n B + l)i {B) p^{t) 

p 34 (t) = -[(n A + l)7 (A) + (^ + ^)7 (i?) ]P34(t)+n A 7 {A) Pi2(t) 

P4i(0 = -[(n A + lh^ + (n B + 1 -h^}p 41 (t) 

p 42 (t) = -[(n A + \h (A) + (n B + l))7 {B) ]p 42 (t) + n B 7 (B) P3i(*) 

p 43 (t) = -[(n A + lh (A) + (n B + ^h {B ^(t)+n A ^p 21 (t) 

pu(t) = -i(nA + lh (A) + (n B + l)7 (B) ]p 44 (t) + n Al {A) p 22 {t) + n Bl ^ B) p 33 (t) 

For the sake of simplicity, we assume that the cavities are identical = 

>y( B ) — ry anc [ fi A — fi B — ji On solving these equations we get the time 

evolution of the density elements matrix 

Pn(i) = ^{[-2 + 2a(2 Pll (0)+p 22 (0) + p33(0)-l)](l + a)e-" t 

+ [l + a 2 (l-2p 33 (0)-2p 22 (0)) 

- 2a(2 Pll (0) + p 22 (0) + p 33 (0) - l)^ 2 "' 

+ (1 + a) 2 } 

p 22 (t) = ^{[-l + a 2 (2p 22 (0) + 2p 33 (0)-l) 

+ 2a(2p n (0) + p 22 (0) + p 33 (0) - l)^ 2 "' 

+ [2 + 2a 2 (p 22 (0) - p 33 (0)) - 2a(2p n (0) + p 22 (0) + p 33 (0) - l))]e"^ 

+ (« 2 -l)} 

P3s(*) = ^{[-l + a 2 (2p 22 (0) + 2p 33 (0)-l) 

+ 2a(2p n (0) + p 22 (0) + p 33 (0) - l)]e^ 

+ [2 + 2a 2 (p 33 (0) - p 22 (0)) - 2a(2p n (0) + p 22 (0) + p 33 (0) - l)] e -* 

+ (« 2 -l)} 

Pu(t) = 1 - pn - p 22 - p 33 

Pi*(t) = ^{[a(pi 2 (0) - p 34 (0)) - p 34 (0) - p 12 (0)]e-i"* 



+ (p 34 (0)+pi 2 (0))(l + a)e-^} 



Pis(t) = ^{[a(-p 24 (0) + pi 3 (0)) - p 24 (0) - p 13 (0)]e-i^ 
+ (p 24 (0)+pi 3 (0))(l + a))e-^} 
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Pu(t) = Pi4(0)e-"* 

P2s(t) = p 2 3(0)e-^ 

Pn(t) = ^{(a(p 2 4(0)-p 1 3(0)) + p 2 4(0) + p 13 (0)) e -| r ' i 

+ (p 2 4(0)+p 13 (0))(-l + a))e-^} 

Pu(t) = ^{[a(P34(0)-p 12 (0))+p34(0) + p 12 (0)] e -|" t 

+ (p34(0)+pi 2 (0))(-l + a)e-^*} 

and p 2 i(t) = p 12 (t),p 31 (t) = p 13 (t),p 32 (t) = p 23 (0,P4i(i) = Pu{t),p A2 {t) = 
P2i(f) 1 Pisit) = p 34 (t), where a = 2n + 1 and r\ = cry. 
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Cavity A Cavity B 




Figure A.l: The scheme for single-particle interference. A two- level atom prepared in its 
excited state passes successively through cavity A, a field damping region and cavity B. 
At the end the two- level atom will be in its ground state. 
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Figure A. 2: Initial value of logarithmic negativity (line) and Concurrence (dot) in terms 
of |a 2 | where |a,2 1 2 + |«3 1 2 = 1 
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;ure A. 3: Entanglement dynamics of the first EPR state for different initial probability 
vacuum reservoirs, (a) |a 2 | = -?=, (b) |a 2 | = 4=. 
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Figure A. 4: Entanglement dynamics of the first EPR state at different mean photon 
numbers of reservoirs in the case |a 2 = |a 3 | = -4=. (a) n = 0.1, (b) n = 0.25. 
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Figure A. 6: Entanglement dynamics of the second EPR state at different mean photon 
numbers of reservoirs in the case |ai| = 04 = -4= (a) n = 0.1, (b) n = 0.25. 
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